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Abstract 

We prove that the asymptotic field of a Skyrme soliton of any degree has a non- 
trivial multipole expansion. It follows that every Skyrme soliton has a well-defined 
leading multipole moment. We derive an expression for the linear interaction energy 
of well-separated Skyrme solitons in terms of their leading multipole moments. This 
expression can always be made negative by suitable rotations of one of the Skyrme 
solitons in space and iso-space. We show that the linear interaction energy dominates 
for large separation if the orders of the Skyrme solitons' multipole moments differ by 
at most two. In that case there are therefore always attractive forces between the 
Skyrme solitons. 



1 Skyrme solitons 

The fundamental field of Skyrme's theory |J is a map 

U : M 3 -> SU{2). (1.1) 

We denote points in M 3 by x with coordinates Xi, i = 1,2,3 and Euclidean length r = 
|x| = + x 2 + x i • Sometimes we write x for the unit vector x/r. It is often useful 
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to parametrise U in terms of the Pauli matrices t% , r 2 and r 3 and the triplet of pion fields 
7Ti , 7r 2 and 7r 2 as 

U(x) = o(x) + i7C a (x)T a , (1.2) 

where summation over the repeated index a is implied and the field a is determined by the 
constraint a 2 + nf + -k\ + 7r| = 1. In this introductory section we do not specify the class of 
functions to which U belongs. It is assumed to be sufficiently smooth for all the following 
operations to make sense. 

The Skyrme energy functional is best written in terms of the Lie-algebra valued currents 

U = U^diU (1.3) 

or 

Ri = UdiU\ (1.4) 

where <9j = d/dxi. It is 

E[U] = -Jd 3 x Qtr(L^) + itra^AHL^Dj . (1.5) 

The Euler-Lagrange equation for stationary points of this functional is conveniently ex- 
pressed in terms of the modified currents 

U = U-\\l 3 A^M\ (1.6) 

and 

Ri = Ri — -g[Rj, [Rj,Ri]], (1.7) 
where we again use the convention that repeated indices are summed over. It reads 

diU = (1.8) 

or, equivalently, 

diRi = 0. (1.9) 

Here we are interested in finite-energy solutions of the Euler-Lagrange equation. It is 
shown in || that the finite energy requirement implies that the map U tends to a constant 
value at infinity in a weak sense. We choose that constant to be the identity element 
1 G SU (2) and demand 

lim U(x) = 1. (1.10) 



The boundary condition (|1.10|) means that the domain of U is effectively compactified 



to a three-sphere. Since the target space is also a three-sphere, maps satisfying (|1.10|) 
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have an associated integer degree. The first rigorous proof that for a finite-energy Skyrme 
configuration in a very general class of functions the degree 



deg[U] = -— J (FxeijktTfaLjLk) {1.11) 



is an integer was given in ||. This result means that the space 

C = {U : M 3 -> SU(2) | E[U] < 00} (1.12) 
of finite-energy configurations is a disjoint union of sectors 

C k = {U : M 3 -> ST/(2) I £[C7] < 00 and deg[U] = k} (1.13) 

labelled by the integers fceZ. 

The symmetry group of Skyrme's theory will play an important role in our discussion. 
The energy functional ( |1.5| ), the boundary condition (|1.10 ) and the degree ( p..ll|) are in- 
variant under the action of the Euclidean group of translations and rotations in M 3 and 
under rotations of the pion fields 

vr a G ab 7r b , G G 50(3), (1.14) 

which we call iso-rotations. Reflections in space S : x 1— > —x and in iso-space n a 1— > — 7r a 
both leave the energy invariant but each changes the sign of the degree. The pull-back of 
Skyrme configurations U via S provides a map 

S:C k ^C_ k , U^UoS (1.15) 

which preserves the energy. 

It was shown in 0] that the energy in each topological sector is bounded below by a 
multiple of the degree. It follows from the results of |J that the bound cannot be attained 
for the standard version of the Skyrme model described here, so that we have the strict 
inequality 

E[U] > 12n 2 \k\. (1.16) 

The bound ensures that the infima 

I k = M{E[U}\U eC k } (1.17) 

are well-defined. The question of whether the infima are attained was first addressed by 
Esteban in the paper @. Amongst other things Esteban proved that, for a suitable class 
of functions, 

h<h + h-i (1.18) 

for all k, I G Z. She also showed that infima are attained provided one assumes the strict 
inequality 

J fc </, + J fe _, (1.19) 
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for all k e Z - {0,±1} and I E Z - {0, fe} in the range \l\ + \k - l\ < \/2\k\. In [[| it 
was shown that the result still holds if one widens the class of allowed functions but the 
inequality ( |1.19| ) remains a necessary assumption in the proof. The strict inequality is also 
of interest in physics. As we shall see it is related to the question of attractive forces in the 
Skyrme model. 

For low values of k the existence and nature of minima of the Skyrme energy functional is 
understood in more detail. For fields of degree one, the highly symmetric hedgehog ansatz 

U H (x) = exp(i/(r) ), (1.20) 

introduced already by Skyrme, leads to an ordinary differential equation for the profile 
function /. With the boundary conditions /(0) = 7r and /(oo) = the resulting Skyrme 
configuration has degree 1 and is called the Skyrmion. It was shown in || that the Skyrmion 
minimises the Skyrme energy functional amongst all degree one configurations of the hedge- 
hog form. In || the existence of the minimum of the Skyrme energy functional in C\ was 
proved, but it has not been established rigorously that minimising configurations have the 
symmetry of the hedgehog field ( |1.20|) . Note that if U\ is a minimal energy configuration 



in C\ then the reflected configuration S (Ui) has the same energy and minimises the energy 
in C-\. 

In the following we use the term Skyrme solitons for minimal energy solutions of the 
Skyrme equation with non-zero degree k. There is overwhelming numerical evidence that 
the minimum in C2 is attained by a configuration of toroidal symmetry 0. For higher 
degree, too, much is known numerically about the minima of Skyrme's energy functional 
in the sector Ck- Numerical searches, assisted by analytical ansatze and investigations of 
the possible symmetries of Skyrme solitons, suggest the existence of Skyrme solitons in all 
sectors Ck up to k = 22. The energies are sufficiently accurately computed that it appears 
that the inequality ( |1.19|) is satisfied for all k in the range 2 < k < 20 and all / in the range 
< I < k, see [§, |, [Lj. 

The existence of attractive forces between two Skyrmions was already shown by Skyrme, 
using the product ansatz. In this paper we investigate the existence of attractive forces 
between general Skyrme solitons. It is perhaps worth stressing that our arguments in fact 
apply to any finite-energy solution of the Skyrme equation, not just minimal ones. An 
earlier attempt at proving the existence of attractive forces between Skyrme solitons was 
made in the unpublished paper [11]. Our approach is partly inspired by ideas in [[nj] but 



also fills important gaps left there. Our main tool is an asymptotic expansion of Skyrme 
solitons. In Section 2 we show that Skyrme fields have an asymptotic expansion in powers 
of 1/r and that for non-trivial Skyrme solitons that expansion always has a non-trivial 
leading multipole. In Section 3 we study the interaction energy of two Skyrme solitons and 
show that it is dominated by a certain linear interaction energy provided the orders of the 
leading multipoles of the Skyrme solitons do not differ by more than 2. In Section 4 we 
derive an expression for the linear interaction energy of two multipoles, and show that it 
can always be made negative by suitable relative rotations in space and iso-space. At the 
end of this paper we briefly comment on the relationship between our results and Esteban's 
work, and on the implications for the existence of Skyrme solitons of arbitrary degree. 



4 



2 The asymptotics of Skyrme solitons 



The aim of this section is to show that if U and the currents Lj are a little better than 
continuous, then U is smooth in M 3 , and has a non-trivial asymptotic expansion in powers 
of 1/r and logr as r — > oo. By non-trivial we mean here that if U is non-constant, then 
it cannot happen that all terms in the asymptotic expansion vanish. Put another way, it 
is not possible for U to approach 1 at infinity faster than every power of 1/r unless U is 
identically equal to 1 in M 3 . 

To make a precise statement, say that the (possibly matrix- valued) function / is in 
Cj' a (R 3 ), with < a < 1, if / satisfies 

mj (l+r + r') a \f(x)-f(x')\ ^ 

SUP f + SUP i : ; r - — < OO. 2.1 

xeR3 x ,x>m( x ^x> \r - r'\ a + (1 + r + r') a d{uj, u') a 

Here x = ru,x' = r'u/, where u is regarded as an angular variable (unit vector) living on 
the unit 2-sphere, and d(u, lj') is the 2-sphere distance between u and to' . If we consider the 
space C^' a {K), where K is a bounded subset of R 3 , then this is precisely the same as the 
usual space of (bounded) Holder-continuous functions, with Holder exponent a. However, 
C°' a (M 3 ) is slightly different from the usual space of bounded Holder-continuous functions 
on R 3 . 

The main result of this section can be summarized as follows: 

Theorem 2.1 Suppose that for some 5 > 0, 

(l + r) 5 (U-l)EC° b ' a , (l + r)%eCj' a . (2.2) 



Suppose further that U satisfies the Skyrme differential equations ( 1.8j) in the sense of 
distributions. 

Then U is smooth in IR 3 and has a complete asymptotic expansion in powers of 1/r and 
logr, for large r. If U is non-constant, then this expansion has a leading term which is 
harmonic, hence a multipole. 



Note that the hypotheses fl2.2|) force U to approach 1 and the Lj to approach like 
r~ 5 as r ^ oo. As a technical remark, we point out that the assumption of Holder- 
continuous currents implies that U will have a Holder continuous first derivative. The 
second derivatives djdkU are then defined only in the sense of distributions, but in the 
Skyrme equation djdkU enters linearly, and is multiplied by continuous functions (smooth 
functions of the currents) see Section |2.1| below. In particular, the left-hand side of the 
Skyrme equation djLj = makes sense as a distribution if ( p.2|) holds. 

The assumptions ( |2.2[ ) of Theorem ( |2.1| ) do not follow immediately from the variational 
analysis used by Esteban in p[. Her methods only give that the derivatives djU are lo- 
cally square-integrable (and that the components of U are locally bounded). On physical 
grounds, one expects minimizers of the Skyrme energy functional to satisfy ( |2.2| ), but it 
would be desirable to bridge the gap between the analysis given here and what was proved 
rigorously in M and B. This issue will not be pursued further here. 



5 



The proof of Theorem ( |2.1| ) proceeds in four steps, each of which takes up one of the 
following subsections. In the first subsection we re-write the Skyrme equation in order to 
make explicit the form of the non-linearities. The equation is quasilinear, in the sense that 
the derivatives of highest order (2) enter linearly. The Skyrme equation can therefore be 
regarded as a second-order linear elliptic PDE, with Holder-continuous coefficients. Then 
standard regularity theorems (Schauder estimates) yield that U is smooth. 

In Section |2.2| we study the behaviour of the Skyrme equation at spatial infinity, by 
introducing coordinates (s,u), where s = r" 1 , so that the 2-sphere at infinity becomes a 
genuine boundary at s — 0. After a rescaling, the Laplacian of M 3 is replaced by 

A b = (sd s ) 2 - sd s + A*,, (2.3) 

where A w is the Laplacian of the unit 2-sphere. The analysis is now guided by the corre- 
sponding analysis of a system of ordinary differential equations with regular singular point 
at s = 0. In particular, one does not expect the solutions to be smooth near s = 0, but 
one does expect a non-trivial expansion in powers of s (and possibly logs). The 6-calculus 



of []T2| enables us to make this kind of argument precise. Thus in Section [T2] we show that 
if U — 1 + u, then u is conormal at s — 0, which is to say that u and all derivatives of the 
form (si9 s ) m V™M are continuous as s — > 0. (This condition is strictly weaker than u being 



smooth near s — 0.) In Section |2.3| , we show that it is not possible for U to approach 1 
faster than r _/i for every fi > unless U — 1 in M. 3 . Finally in Section we combine 
this fact with another application of the 6-calculus to show the existence of a non-trivial 
asymptotic expansion in powers of r~ x and logr (or equivalently in powers of s and logs). 

2.1 Rewriting the Skyrme Equation 

To begin with we work near a fixed point of M 3 , which we may as well take to be the origin 
0. By replacing U(x) by f/(0) _1 f/(x) we can assume that U(0) = 1. Write 

U(x) = l + u(x), (2.4) 

so that u(0) = 0. Because U is unitary, the 2x2 complex matrix u will satisfy the algebraic 
constraints 

u + u ] + uu ] = tr(tt) + detw = 0. (2.5) 
In particular u is neither exactly skew-hermitian nor trace- free. Then 

Li = diU + u^diU (2.6) 

and 

diU = (1 + u ] ) Am + diridiU. (2.7) 
The cubic term in the currents can be written 

-Lj[Lj, = ^{dju + vtdjii)[djU + u^dju, diU + v)diu}. (2.8) 

Taking the divergence, and using the notation 

Vij = diU^djii, (2.9) 
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we obtain 



where 



l -L ] [L J ,L l } = {l + u^){T + F), (2.10) 



and 



T = T{u, du, d 2 u) = _tH (1 + u + )Au] + Lj[{\ + v))diud jU , L t ]) (2.11) 

1 + u 

F = F(u, du) = —^-( Vij [Li, Lj] + Lj[Lj, v u ) + L 3 [(l + v))v ij: U}). (2.12) 

The nonlinear terms have been divided here so that T is polynomial in the first derivatives 
of u and linear in its second derivatives, while F only contains u and its first derivatives. 

With this notation, the full Skyrme equation can be written 

P(u, du, d 2 u) = Q(u, du) + F(u, du) (2.13) 

where 

P(u, du, d 2 u) = Au + T{u, du, d 2 u) and Q(u, du) = -(1 + u)v u . (2.14) 

Notice that Q is (approximately) quadratic in u, T is of degree three and F is of degree 
four. Note also that T is linear in didjU (and quadratic in d{u). 

Because of the assumed Holder continuity of the currents, the coefficients of the differ- 
ential operator 

f^P(u,du,d 2 f) (2.15) 

are Holder continuous, and this operator is linear and elliptic in a small neighbourhood K 
of 0. Moreover, the RHS Q + F of ( 2.13|) is also in C®' a (K), so by elliptic regularity, it 
follows that u G C b ' a (K). Then the currents are in C^ a {K) and the process continues to 
show that u is in C^' a (K) for every k. Thus u is smooth near 0. Since the point was 
arbitrary, the argument shows that u is smooth in M 3 . 



2.2 Boundary regularity 

In order to analyze the asymptotic behaviour of the field U, we shall make a transformation 
of the problem which involves passing from M 3 to a 'compactification' R 3 in which the 
sphere at infinity becomes a genuine boundary. This is easily achieved by introducing the 
coordinate s = 1/r along with angular coordinates 9 and tp in M 3 . Then in the space 
[0, oo) x S 2 with coordinates (s, 9, <p), the boundary s = corresponds to r = oo and 9 and 
(p give coordinates on the boundary, which is the 2-sphere 'at infinity'. M 3 is defined to be 
the union of M 3 with the 2-sphere at infinity attached in this way. 

It can be cumbersome to work with explicit coordinates on the 2-sphere, so we again use 
uj for points on S* 2 and represent any point other than the origin of M 3 in the form (s, uj). 

Next, introduce rescaled derivatives 

Di = Tdi=-di. (2.16) 
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These vector fields have the property that they are linear combinations, with coefficients 
that are smooth, down to s = 0, of the basic vector fields sd s , do and d^. The Euclidean 
Laplacian takes the form 

A = s 2 A b (2.17) 



with A b defined in (|2.3|). 

Now we write U = 1 + u for large r (that is, for small positive s) and make the rescalings 
of ( |2.13| ) suggested by ( [2.16[ ) and ( |2.17| ). The result is a '6' version of the Skyrme equation, 

P b {u, Du, D 2 u) = Q b (u, Du) + s 2 F b (u, Du), (2.18) 

where 

P b (u,Du,D 2 u) = A b u + s 2 T b (u,Du,D 2 u), (2.19) 

Q b (u,Du) = -(l+u)D iU j DiU, (2.20) 

these terms being got by replacing di by the rescaled derivative Di wherever they occur. 

The reason for reformulating the equation in this way is that there is a well-established 
theory, called the 6-calculus, which can be used to analyze equations of this kind |L2| . The 
6-calculus is concerned with 6-differential operators. In the present situation, a ^-differential 
operator is just a differential operator of the form 



P= J2 C abc {s,u){sd s ) a d b e d^ (2.21) 



a+b+c<m 

where the coefficients C abc are smooth up to the boundary s = 0. From ( |2.3|) it is clear that 
the rescaled Laplacian A b is an example of such an operator. The set of all such operators 
will be denoted by Diff fe , those of order k by Diff^. 

The first aspect of this theory that is needed is the counterpart of the elliptic regularity 
for Holder spaces that we used in the previous subsection. For this we need 6-H61der spaces, 
already introduced in fl2.1|) . With the new variable s, we have f{s,u) e C°' a if 

(s + s') a \f(s,oo)- f(s',u')\ 

s »p / + su p r — /ill \ Jm fil <0 °- 2 - 22 

Now put 

a = {f : Lf E C° b ' a for all 6-differential operators L G Diff£} 

and 

A = {u : Lu E C° b ' a for all L G Diff b }. 

In order to force functions to decay as s — > 0, we introduce weighted versions of these 
spaces, 

s 5 C^' a = {u = s 5 v:vE Ct'"}, s 5 A = {u = s 5 v : v G A}. 



The following is a very special case of elliptic regularity for 6-differential operators [13 
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Theorem 2.2 Consider the differential operator P = + s 2 E, where E is a second- 
order differential operator with coefficients smooth up to the boundary s = 0. Suppose that 
Pu = f near s = 0, with u and f in s s C®' a . Then if 5 is not an integer, it follows that 
ues 5 C 2 b a . 

We want to apply this to the Skyrme equation, written in the form ( |2.19[ ). After the 
last section, we know that the coefficients of the perturbing term E are smooth for s > 0, 
but all we know at s = is the original assumption that the currents are in s s C®' a . The 
elliptic regularity statement still holds in this case (though this statement does not seem 
to be available in the literature). Applying this result, for 8 a small positive number, we 
obtain first that u G s 5 C 2,a , which gives that the coefficients are in s s C b ' a , so u G s s C b ,a 
and so on. Iterating, we find u G s s A. 

This is a major step forward: in particular, the angular variation of u is now very well 
controlled. However u may still be very far from being smooth up to the boundary or 
having an asymptotic expansion there. Indeed, horrors like s s sin logs lie in s s A. 



2.3 Power-law decay of u 



In this section and the next we establish that a topologically non-trivial solution of the 
Skyrme equation must have a non-trivial asymptotic expansion in powers of 1/r. We show 
first that it is not possible for a topologically non-trivial solution to approach 1 faster than 
every power of 1/r. This result is then fed into an iterative analysis of the equation in the 
next section. These two sections can, however, be read in either order. The main result of 
this section is as follows. 

Theorem 2.3 If the topological charge of the Skyrme field U is non-zero and if it satisfies 
the hypotheses of Theorem \2. 1\ then there exists some fi G M + such that r M w(r, uj) does not 
tend to zero as r — > oo. 

To set this result in a more general context, recall that a partial differential equation 
Lu = is said to have the unique continuation property at a point 0, say, if the following 
is true: 

If all derivatives of u vanish at 0, then u = in some neighbourhood of 0. (2.23) 

The methods used to establish that a second-order PDE has the unique continuation prop- 
erty establish analogous statements with somewhat weaker hypotheses. For example, a 
simplified version of Theorem 17.2.6 of |1| is as follows. 



Theorem 2.4 Let ajk(x) be smooth and positive-definite in a neighbourhood X of 0, and 
suppose that 0^(0) = 8^. Suppose that for x G X, the smooth function u satisfies 

{ajkWdjdku] < A (\u(x)\ + \Vu(x)\) (2.24) 

for some constant A, and 

|x| _m |m(x)| — >• as \x\ — > for every /i G M + . (2.25) 
Then u = identically in a neighbourhood ofO. 
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Using this result we can prove the following: 



Theorem 2.5 Let U satisfy the hypotheses of Theorem ( \2.1\ ). Then if the topological charge 
of U is non-zero, U cannot be constant in any open subset ofM. 3 . 

Proof Define 

W = {x G M 3 : U(y) = U{x) for all y in some neighbourhood of x}. (2.26) 
Then W is open by definition. By Theorem [2.4|, W is also closed. To see this, suppose that 



x n G W, x n — > Xq G M 3 . Then all derivatives of U are zero at Xq, so ( |2.25| ) holds (with 
replaced by xq). The Skyrme equation, in the form ( |2.13| ), implies a differential inequality 
of the form ( |2.24|) in a neighbourhood of xq- It follows^] that U is identically constant in a 



neighbourhood of x , so that x G W. Since M 3 is connected, W = or W = M 3 . □ 

We will use this theorem to give an indirect proof of Theorem |2.3| . Suppose r^u(x) — > 
as r — > oo for every /x. Adapating Theorem p.4| we shall show that then u = for 
a// sufficiently large r. Thus {7 is constant in an open set, hence by the previous result 
constant everywhere. 

So consider the 6-differential operator 

P b = A b + s 2 T b (2.27) 



on R 3 , where the coefficients of T b are in C b ' a . By rescaling Theorem 17.2.6 in |16 
obtain 

Theorem 2.6 Suppose that 

\P b u(s,cj)\ < As 5 {\u{s,uj)\ + \Du(s,u)\) for all < s < s ,u G S 2 (2.28) 

and that 

s""|u(s,a;)| -> as s -> /or a// //. (2.29) 
T/ien tt(s, cj) = /or < s < si, where si is some small positive number. 



The '6' version ( p. 18 ) of the Skyrme equation implies a differential inequality of the form 



(|2.28|) , just as before. It follows that if u decays faster than any power of r, then u = for 



all sufficiently large r. By the remarks before the statement of Theorem |2.6| , the proof of 
Theorem |2.3| is now complete. 



1 Unique continuation theorems do not extend wholesale to systems. However, in our case, the leading 
term is the Laplacian and the other second-order terms C(u,du,d 2 u) are non-scalar but small near xq. 



The proof of Theorem 2.4 goes through in this case. 



10 



2.4 Refined regularity, asymptotic expansions 

The main result of this section can be stated as follows: 



Theorem 2.7 Let U = 1 + u satisfy the hypotheses of Theorem \2. 1 . Then there is some 
integer M > 1 and an asymptotic expansion 

2M oo 

u ~ Yj(Ld)r~ti +1 ^ + r~^ +1 ^Wj(u, logr) for large r. (2.30) 

j=M jr'=2A/+l 

Here the Yj are Lie-algebra valued spherical harmonics, 

X)) ,/U • uv,, (2.3i) 



and Ym ^ ; so that the piece Y^j=m ^j ( Co ') r a non-zero harmonic function. The 

functions Wj are smooth in uj and polynomial in logr. 

It will follow from the proof that the terms u?2A/+i to w^m+i are of at most degree 1 in 
logr, the terms wsm+2 to are of at most degree 2 in logr and so on. We remark also 

that the asymptotic expansion can safely be differentiated term by term to give asymptotic 
expansions of all derivatives of u. 

In order to motivate the proof, consider the equation 

A b u = f (2.32) 

where u and / are defined for small s. If / has the form / = s x g(u), then a solution can 
be found as follows. Expand g as a sum of spherical harmonics, 

oo 

9 = ^29j, where A^gj = + l)gj 

j=0 

and seek a solution u = Y2 u j( s )9j- Then Uj must satisfy 

[(sd s ) 2 -(sd s )-j(j + l)} Uj (s) = s\ 

This is solved by 

provided there is no resonance, that is to say 

X^-j,j + l. 

In the resonant case, Uj(s) has the form s x (A + B logs). The general solution is obtained 
by combining this with an arbitrary solution of the homogeneous equation A^v = 0. If we 
require u — > as s — > 0, then v must itself go to zero and hence will be a sum of multipoles, 
v = YlJLo^ji 00 ) 3 ^ 1 ' wnere Yj satisfies Q2.31|) . 
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The 6-calculus extends results of this kind to functions in s s A (which, as we have seen, 
can be far from having expansions in powers of s). In order to summarize the needed 
results, write 

/ = 0(s 5 ) instead of / G s 5 A. (2.33) 

This conforms to the use of the 'O'-notation in the rest of the paper, but has the additional 
property 

If / = 0(s 5 ) then Lf = 0(s 5 ) for all L G Diff 6 . (2.34) 



Lemma 2.8 Suppose that u and f defined near s = satisfy dPID - Ifu = 0(s 5 ) and 
f = 0(s n+s ), where 5 > and n is a positive integer, then 

n-l 

u = Y i(u)s j+1 + 0{s n+s ) (2.35) 

3=0 



where Yj satisfies ($.3\ ). In particular the sum on the RHS is harmonic 



n-l 

A b J2Yj(u)s j+1 = 0. (2.36) 

3=0 



Lemma 2.9 Suppose that u and f defined near s = satisfy ft2M ). If u = 0{s s ) and 



n-l 

f = ^/>,log S y +1 + 0(s n+5 ), (2.37) 

3=0 

where 5 > 0, n is a positive integer and fj is polynomial of degree rrij in logs, then 

n-l 

u = ^u^,logsy' +1 + 0(s n+s ) (2.38) 

3=0 

where Wj is a polynomial of degree rrij + 1 in logs. 

These results will be applied to the Skyrme equation, now rewritten as 

A b u = Z b {u) := Q b {u) - s 2 T b (u) + s 2 F b (u). (2.39) 

Lemma 2.10 Suppose that U = 1 + u satisfies the Skyrme equation and u = 0(s s ). Then 
u = 0{s 2 ). 

Proof: On the RHS of (|2.39|) , Q b is quadratic in Du and the other terms in Z b are of even 
higher degree. Hence Z b = 0(s 25 ). Applying Lemma |278| we obtain that u = 0(s 2S ) + 0(s). 
If 25 < 1, we can iterate this argument to obtain eventually u G O(s). In ]14|]> it is 
shown that a solution of the Skyrme equation cannot have leading term 1/r = s in its 
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asymptotic expansion. It follows that u = 0(s 1+s ) (for a possibly smaller 5 > 0). Hence 
Z b {u) = 0(s 2+2S ) and so, applying Lemma |2]8| again, u = Yis 2 + 0(s 2+s ). □ 

Combining Theorem |2.3| with Lemma ET8L we see that there exists an integer M > 1 with 
the property that 

u = Y M {u)s M+1 + 0(s M+1+8 ), (2.40) 

where Ym is a non-vanishing spherical harmonic. We can now complete the proof of Theo- 
rem 2]7 in the following iterative fashion. From the structure of Z b (u) it follows from (|2.40|) 
that 

Z h {u) = f 2M+1 s 2M+2 + 0(s 2M+2+s ). (2.41) 



Applying Lemma 2.9, 



u 



2M 

YA"> j+1 + w 2M+1 (u,logs)s 2M+2 + 0(s 2M+2+s ), (2.42) 

j=M 



where W2M+1 (^, log s) is of degree at most 1 in log s. Now this expression for u is substituted 
into Z b) to give 



3A/+1 



Z b = J2 fj("> j+1 + f3M + 2(oo,\ogs)s 3M+3 + + 0(s 3M+3+5 ) (2.43) 



2M+1 



where /3M+2 is of degree at most 1 in log s. Now apply Lemma |2l] to get 

2M 3M+1 

u= ^7 J ( W y +1 + ^^gs)s' +1 + w 3M+2 (u,logs)s 3M+3 + 0(s 3M+3+5 ),(2AA) 

j=M j=2M+l 

where the functions W2M+1, w ?j m+i are of degree at most 1 in log s and W3M+2 is of degree 
at most 2 in log s. Carrying on in this way we obtain a complete asymptotic expansion. 

To complete the proof of Theorem p.7| note finally that the first M + 1 terms in the 
expansion are actually the pion field. Indeed, if the expansion ( |2.30| ) is substituted into the 
algebraic constraints (|2.5|) we see that the harmonic piece J2"m ^j' 5 ^ 1 ^ s skew-adjoint and 
trace- free — the quadratic corrections enter at order s 2M+2 . □ 

The upshot of this section is that every Skyrme soliton has a leading Lie-algebra valued 
multipole field (called a 2 M -pole) 

M 

»»W= ir -E^I«.^ (2-45) 



where Yum are the usual spherical harmonics on S 2 , see appendix A. The leading multipole 
moment Q a Mm is independent of the location of the Skyrme soliton, and is acted on naturally 
by rotations and iso-rotations. It is a key ingredient in the calculations of the following 
sections. As already mentioned one can show that Skyrme solitons never have asymptotic 
monopole fields ||14|| . The leading multipole field of the B = 1 hedgehog (|1.20| ) is a triplet 
of dipoles, and dipoles are known to occur as leading multipoles in a number of Skyrme 
solitons. The highest leading multipole known from numerical work is an octupole, which 
occurs in a B = 7 configuration with icosahedral symmetry B. 
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3 The interaction energy of two Skyrme solitons 



Suppose we have Skyrme solitons and of degrees k and I which minimise the 
energy in the sectors Ck and Q. Since the total energies of and are finite there 
must be balls B\ and B 2 in R 3 so that most of the energy of and is concentrated 
in, respectively, Bi and B 2 . Outside the balls B\ and B 2 the asymptotic analysis of the 
previous section applies. Suppose that the leading multipole of is a 2 M -pole and the 
leading multipole of is a 2 Ar -pole. Denoting the radii of Bi and B 2 by D\ and D 2 , and 
with the abbreviation (|2.45 ) for a generic Lie-algebra valued multipole field we have 



U {1) (x) ~ 1 + u M (x) for x <£ B x (3.1) 

and 

U [2 \x) ~ 1 + ^(x) for x £ 5 2 . (3.2) 

Using the translational invariance of the Skyrme energy functional we can assume without 
loss of generality that B\ is centred at X + = (0, 0, R/2) and that B 2 is centred at X_ = 
(0,0, —R/2), where R is so large that B\ and B 2 do not overlap, i.e. R > D\ + _D 2 . The 
parameter R will be interpreted as the separation of the Skyrme solitons. There is clearly 
an ambiguity in the definition of such a separation parameter, but this does not affect 
our calculation of leading terms in the limit where R becomes large. Then we define the 
following product configuration 

U R (x) = U {1) {x)U {2 \x). (3.3) 

This configuration has degree k + / and we shall see shortly that its energy is finite, so that 
Ur £ Ck+i- 

Our goal is to study the energy of the product configuration Ur perturbatively in the 
limit of large R and to compute the leading terms in powers oil/ R. A similar calculation 
for moving and spinning Skyrmions was performed in ||15|| , where some further details are 
given. Let uf 3 and Lf' be the currents ( |1.3|) and (|1.6| ) constructed out of and R^ 
and R^ be the currents ( p..4|) and (\1.7j ) constructed out of . Then one computes 

E[U R ) = E[U {1) ] + E[U {2) ] + W 2 + W 4 . (3.4) 

The energies EpW] and E[U {2) ] are simply the energies of the Skyrme solitons U^ 1 ' and 
and therefore independent of R. The interaction terms W 2 and W& are given by 
integrals over M 3 

W 2 = J d 3 xw 2 and W4 = J d 3 xw 4 , (3.5) 

with integrands 



«*=te(i,«fl« + W-L{ 1 >Jjn (3-6) 
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and 

«. = -|tr ([Lf», fif >][L« flf »] + [Lf», flf][flf», if']) + [Lf», L«][flP>, flf ]) . (3.7) 

We shall see shortly that the term W 2 contains the leading contribution to the interaction 
energy. However, for the presentation of our method of computation it is more convenient 
to begin with the quartic term W4. We split the integration region R 3 into the balls Bi 
and B 2 and the complement C = M 3 — (Bi U B 2 ). To illustrate our method, consider the 
integral 

I = -\ ^ d 3 ,tr ([L?\l$W,Bto]) < \ cfer tr (flf if ) . (3.8) 

The currents Lp are smooth functions and hence bounded on the compact domain B\. 
Therefore 

|i| <~K J d 3 xti (RfRf^j . (3.9) 

for some positive constant K. Since -Bi is far away from the centre of soliton , the lead- 
ing contribution to the integral ( |3.9| ) is obtained by replacing Rj by the leading multipole 
component —idjVN- Using 

for a further positive constant K' we conclude that 



A similar calculation for the other terms in W4 shows that 



B 



*™* = °[j$N+l)- ( 3 - 12 ) 



Considering the contribution from B 2 we find by the same argument 

j B f xw ^ = °{w^)' (3 ' 13) 

Finally the remaining integral over C can be estimated as follows. Define 

nR) = J/* lx _l r ^ x _l r+ „ (3.14) 

noting that the integral converges for all values of R > D\ + D 2 . Then there is a positive 
constant K" such that 

/ d 3 xw 4 < K" F(R). (3.15) 
Jc 
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The large R behaviour of F can be estimated by a scaling analysis. Fix R > D\ + D 2 and 
consider R > Rq. Changing integration variables x i— > (Rq/R)x one computes 



F(R) = ( 



D \ 2M+2jV+5 
Kq 



FGRoH / d 3 x ' 1 



i i~ y I2M+4 l~ y I2V+4 

, S 1 (R)US2(R) F — A +l |£ — A_| 

where Si(R) and S 2 (R) are the thick shells 



(3.16) 



5 2 (i2) = [iGtfKV^^ |(xi,X2,x 3 + CRo/2)| <^ 2 } (3.17) 
which converge to punctured balls 

B? = {x 6 1R 3 - {(0, 0, -Ro/2)}| x 2 , x 3 - {Ro/2)\ < D{\ 

B° 2 = {xeR 3 -{(0,0,-R /2)}\\(x h x 2 ,x 3 + (R o /2)\<D 2 } (3.18) 

in the limit R — > oo. In that limit the integral over Si(R) diverges like R 2N+1 and that 
over S 2 (R) like R 2M+1 . Combining this with with the factor ^-(2M+2Ar+5) we ^ ec [ uce that 
F(R) decays for large R like R~( 2N + 4 ) and ir (2M+4) , just like the contributions ( gig ) and 
(|3.13|) . Combining all the terms, we conclude that the leading terms in W4 decay for large 
R according to R'^ 2N+ ^ and _R-( 2M + 4 ). 

In order to evaluate W 2 we first divide the region of integration into the half-spaces 

H + = {x e R 3 I s 3 > 0} and H = {x e R 3 | x 3 < 0}. (3.19) 

In H + we replace by the leading term and in H~ we replace by the leading 

contribution 1 + um- The result is 

W 2 w ^^itr^^J^Ar,^]]) -tr^^) 

- J d 3 x itr (iiS^i^, W]]) + tr (Wfl? ) . (3.20) 

Now integrating by parts and using the Euler-Lagrange equations diL\ = diR? = for 
the individual Skyrme solitons we convert two of the terms into an area integral 

-/ d 3 xtx (diVNL^ + I d 3 x tx (diU M Rf ] ) 

= / dx\dx 2 ti (vnL^ + umR^} ■ (3-21) 

Since the xiX2-plane is far away from both Skyrme solitons the leading contribution to this 
area integral can be expressed entirely in terms of the asymptotic fields: 

3 

AE = 2 / dxidx2 (um&Vn ~ v%dau a M ) . (3.22) 

0=1 ^3=0 
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A simple scaling analysis shows that AE falls off like R-( N + M + 1 ) f or large R. We skip the 
details here because we shall show how to evaluate AE exactly in the next section. The 
remaining terms in ( |3.20| ) can be estimated with the techniques used in estimating W^. 



The result is 



W 2 = AE + O + O (j^—) (3.23) 



Combining all terms in (|3.4j) we conclude that 



E[U R ) = E[UU] + EpW] +AE + O + O . (3.24) 

Note that AE is the leading contribution to the interaction energy if \N — M\ < 2 i.e. if 
the orders of the leading multipoles of the two Skyrme solitons differ by at most two. We 
will comment on the validity of this assumption at the end of this paper. 



4 Harmonic functions and their interaction energy 

In order to compute the interaction energy AE we need to derive some general results 
about harmonic functions. We define the regions 

Hg = {x e R 3 | x 3 < 5} and Hj = {x e R 3 1 x 3 > -5}, (4.1) 

where the positive parameter 5 is introduced for technical reasons. Then we introduce the 
spaces 

H- = {f:Hi^R\Af = 0, lim/(x) = 0} (4.2) 

r— >oo 

and 

H + = {g : H+ -»• E | A# = 0, lim g(x) = 0}. (4.3) 

r— >oo 

Elements of 7i~ tend to zero at the boundary "at infinity" of , elements of 7i + tend to 
zero at the boundary "at infinity" of Hf. No additional restriction is placed on behaviour 
at the boundaries x 3 = ±5. 

For the calculations in this section it is convenient to split M 3 into K 2 xl and denote 
vectors in R 2 by bold letters, e.g. x = (xi,X2)- We then write x = (x,x 3 ). The most 
general element of Ti~ can be written as 

f d 2 k 

f(x) = J J7^y^ P( k ) exp(ik-x + kx 3 ), (4.4) 

where k = Vis 2 ' and the volume element d 2 k/ ((2it) 2 2k) arises from the combination of d 3 k 
with the delta-function 5(k 2 — k 2 ) which ensures that / satisfies the Laplace equation. Since 
/ is real the Fourier transform p satisfies 

p(k)=p{-k). (4.5) 
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Similarly, the most general element of 7i + can be written as 

f d 2 l 

with I = Vl? and 

q(k)=q(-k). (4.7) 
There is a natural pairing between elements of the 7i~ and those of 7i + 

(/, 9) = [ d Xl dx 2 {gd 3 f - fd 3 g) for feH^gE H + , (4.8) 

Jx 3 =0 



Using the expansion ( |4.4j) and (O) we find, in terms of the Fourier modes, 



This pairing is of interest to us since the interaction energy fl3.22|) is proportional to the 



sum over the pairings (u%j, v%). Therefore, we also refer to the expressions ( |4.8| ) and ( |Q|) 
as the interaction energy of the harmonic functions / and g. 

It is clear that the pairing (|4.8| ) may vanish for some pairs of harmonic functions / G 
g G 7i + . This happens for example if the support of the Fourier transform p is 
complementary to that of the Fourier transform q. However, we shall now show that 
the interaction energy of multipoles can always be made non-zero by rotating one of the 
functions. 

4.1 Multipole fields 

For / G Hr we have an alternative expansion in spherical harmonics Yum 
ft \ 4?r n Y Mm (0 + ,(p + ) 

/(*) = !, 2^ f2M+TT w — ^ — ' ( } 

where r + = \x — X + \ and (9 + , ip + ) are spherical coordinates centred at X + = (0, 0, R/2). In 
this section we only need to assume R > 0, but in our applications we will be interested in 
the large R limit. The coefficients Qum, are called the multipole moments of the function 
/. Assume that / has non-vanishing multipole moments and suppose M is the smallest 
integer such that Qum ^ for some m = —M, M. The function 

/m(x)= (2M + 1) ^ QMm ^ • 

V ^ I m=-M + 

is a 2 M -pole field and Qum are the leading multipole moments of /. 
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It is often convenient to write multipole fields in terms of partial derivatives of the 
Coulomb potential centred at X + : 

(f) + (x) = — = 1 , (4.12) 

r+ y/fP + (x 3 -R/2)*' 

where p 2 = x 2 . The function d^dpd? 1 <j>+(x) is a 2 M -pole field if m 1 + m 2 + m 3 = M. 
However, not of all of the fields obtained in this way are independent. We introduce the 
complex derivatives 

d=hd 1 -id 2 ) and d=^(d 1 +id 2 ) (4.13) 
and note that A = d 2 + Add. Then, since A0 + (x) = we have 

ddMx) = -\d 2 3 Mx)- (4-14) 

Thus a basis for 2 M -pole fields is given by d™ 3 d n 8 n 4> + (x) where M = m 3 + n + n and either 
n or n can be taken to be zero. In appendix A we derive the exact relation between the 
functions <9f +m d~ m (/) + (x), -M < m < 0, and <9f - m <9 m 0+(x), < m < M, on the one 
hand and the spherical harmonics Yum centred at X + on the other. The result is that we 
have the alternative expansion of an 2 M -pole field 

f M (x)= Yl A Mm dl [+m d- m ( t )+ {x)+ J2 A Mm d 3 M - m B m (t )+ {x) ) (4.15) 

-M<m<0 l<m<M 

where the coefficients A Mm , m = —M,...,M are directly proportional to the multipole 
moments Qum- It follows from the results in appendix A that 



Air " (_-\\M+mom 

--Qum (4.16) 



2M+ 1 y/(M -m)\(M + m)\ 



for m > and 



, 4tt (-l)M 2 IH 
2M + 1 yJ(M - m)\(M + m)\ 

for m < 0. Note in particular that the reality of Jm is equivalent to 

A-M(-m) = A Mm . (4-18) 

Multipole fields have a remarkably simple Fourier transform, which will be important for 
us. We use the representation 



x 



/d 2 k r 
— e -k\x3-s\ exp (ifc- x), (4.19) 
2ixk 
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which can be verified as follows. Exploiting the invariance of <p + under rotations in the 
XiX2-plane we may assume that x = (p, 0). Using polar coordinates (k,i[>) for k we first 
carry out the dk integration and then the angular integration: 

27f# f °°dke- k ^-^ex V (ik-x) = I' '^Z- — 1 7T77Ti\ ( 4 -20) 

,(4.21) 



2ti J 1 J 2tc (ipcosip — \xs — R/2\) 

dw -2 



s i 2ni (ipw 2 — 2|x 3 — R/2\w + ip) 

where we changed variables to w = e 1 ^ in the last line. Expanding the integrand in partial 
fractions and using the residue theorem then yields the expression ( |4.12j ). We compute the 
Fourier transform of the multipole field (|4.15| ) by differentiating ( }4.19|) under the integral 
sign. Note that 

d exp(ik ■ x) = -ke~ 1 ^ exp{ik ■ x) and d exp(ik ■ x) = -ke 1 ^ exp(ik ■ x), (4.22) 
so we find 

f M (x) = — / d 2 k e~ k ^~^ k M ~ l U Am*.*"* exp(zk ■ x). (4.23) 

71 J -M<m<M ^ ' 

Here we have used that £3 < S so that in particular x 3 < y. Thus with the normalisation 
we arrive at the following simple expression for the Fourier transform 



(• \ \ m \ 
\) A Mm e im t (4.24) 

Remarkably, this function factorises into a /c-dependent part and the function 

(• \ \ m \ 
\ ) A Mm e im ^ (4.25) 

of the angle ■0. The fc-dependent part e~ k ^k M is non-zero for k 7^ and the function 
only vanishes identically if Aum = for all m = —M,...,M, i.e. if the 2 M -pole field is 
trivial. 

4.2 The interaction energy of two scalar multipoles 

The interaction energy of two multipoles can be expressed in a remarkably compact way. 
Let 

4>-{x) = (4.26) 

be the Coulomb potential centred at X_ = (0, 0, —R/2) and consider the multipole field 
9n{x)= J2 B Nn d^ n d~ n ^{x)+ B Nn d^ n d n ^{x) (4.27) 

-7V<n<0 Kn<N 
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with Sjv(-n) — Bx n . By the same calculation as for Jm above we find the Fourier transform 
of in the XiX2-plane: 

(• \ M 
~) B Nn e m ^. (4.28) 



The interaction energy of the two multipole fields /m and gN 

V MN =(f M ,g N ) (4.29) 



can now be computed using the formula ( |4.9| ). Using the factorisation property of the 
Fourier transforms Pn and qw, it is easy to perform the integration over k. Assuming 
without loss of generality that M < N we first carry out the integration over the angle ijj 
to find 

poo M 

Vmn = 4vr / dke~ kR k N+M ^ 2~ 2 M A Mm B Nm , (4.30) 

J ° m=-M 

where we have used the reality condition for the coefficients Aum, and B^ m . Computing 
the remaining integral we obtain the final result 



vmn = Air m+n+1 } j 2 1 l A Mm B Nm . (4.31) 



M 



ftM+N+1 



This formula has a number of interesting features. The interaction energy depends only 
on the separation of the multipoles and on the combination J2m=-M 2~ 2 ^A Mm B Nm of 
the multipole components. As explained in appendix A, the multipole moments Qnu of 
a 2 Ar -pole can be thought of as elements of the (2iV + l)-dimensional unitary irreducible 
representation Wn of 5*0(3). The vector B with 2N + 1 components B^n, ~N < n < N 
is naturally an element of Wn- Rotations G G 50(3) about the centre A_ of the multipole 
field qn act on the multipole components via B^n i— > Yln'=-N ^nn'(^)BNn', where U N is a 
(2iV + l)-dimensional irreducible representation of 50(3) (because of the rescaling (|4.16|) 
and (14. 17 ) this is not the standard unitary representation). With our assumption that 



M < N we can use the multipole components AM m , —M < m < M, to define the linear 
form 



M 



F A :W N ^R, B^ 2- 2lml A Mm B Nm . (4.32) 



-M 



By assumption, the components Au-m are not all zero, and therefore the map F A is non- 
degenerate. Writing the formula (|4.31|) in terms of this map as 



(M + N)\ 

Vmn = ^ K rM+n+ [ F A (B), (4.33) 



we immediately deduce the following result. 
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Theorem 4.1 The interaction energy of an 2 M -pole and a 2 N -pole separated by a distance 
R is always non-vanishing for some relative orientation of the two multipoles. When such 
an orientation is chosen, the modulus of the interaction energy decreases with the separation 
asR-( M + N+1 l 

Proof: Assuming without loss of generality that the multipoles are separated along the 
a;3-axis and that M < N we have the formula ( |4.33| ) for the interaction energy. Since the 
map ( 4.32 ) defined in terms of the (non-vanishing) multipole components Aj^ m of the 2 M 



pole at X + is non-degenerate it has a 2iV-dimensional kernel. It then follows from the 
irreducibility of the (2iV + l)-dimensional representation Wn that there exists a rotation 
G G SO (3) such that U N (G)B is not in the kernel of Fa for some G. For that G we thus 
have F A (U N (G)B) = « ^ and V MN = 4ttk(M + N)\R^ M+N+1 \ □ 



5 Attractive forces and existence of minima 

The arguments of the previous section apply to the asymptotic pion fields of the Skyrme 
solitons and discussed in Section 3. In particular we note that the interaction 
energy AE ( |3.22j ) for the leading multipole fields um an d is just the sum over iso- 
components of pairings of the form ( |4.29| ) 

AE a = -2(u a M ,v%) = 2 I dx x dx 2 (u a M d 3 v% - v a N d 3 u a M ). (5.1) 

Jx 3 =0 

Now pick one of the iso- indices, say a — 1, and use iso- rotations of the Skyrme solitons to 
make sure that the first iso-components v} M and v x N are non-vanishing. It then follows from 
theorem ( |4.1|) that we can make the multipole interaction energy AEi non-zero by spatial 
rotations of one of the Skyrme solitons. This fact is the main upshot of the calculations in 
the previous section and a crucial input for the following argument which was missing in 

Now consider the sum 

AE = AE 1 + AE 2 + AE 3 . (5.2) 

We would like to show that we can always arrange for AE to be negative by a suitable 
iso-rotation of one of the Skyrme solitons. We may assume that, possibly after re-labelling 
the pion fields, 

AEi > A£ 2 > A£ 3 , (5.3) 

If AE < we are done, so suppose that AE > 0. Since we know that not all AE a vanish 
we can conclude that AE\ > 0. Now perform an iso-rotation of Skyrme soliton 1 by 180 
degrees around the third iso-spin axis. This reverses the sign of 7r[ and tt^ and hence 
also the sign of AEi and AE 2 . The new value of AE is 

AE = —AEi - AE 2 + AE 3 

= -AEi - (AE 2 - AE 3 ) < (5.4) 
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since —AE\ < and, with our ordering, — (AE 2 — AE 3 ) < 0. 

Thus, the contribution AE to the interaction energy of two Skyrme solitons and 
can always be made less than zero by suitable rotations and iso-rotations of Skyrme 
soliton 1. It follows from the discussion at the end of Section 3 that for \N — M\ < 2 and 
sufficiently large separation parameter R, 

E[U R ]<E[U^} + E{U {1) }. (5.5) 

We conclude with a few comments on the implications of our result for the question 
of existence of general Skyrme solitons. As explained in Section 1, Esteban proved the 
existence of Skyrme solitons of arbitrary degree provided the strict inequality ( |1.19| ) holds. 
Our result ( |5.5| ) implies the inequality in those cases where minima exist in the sectors I 
and k — I, and where the associated multipoles have orders which do not differ by more 
than two. Since monopole fields do not arise in Skyrme solitons, the interaction energy 
AE dominates at large separation if the leading multipole moments in Skyrme solitons 
are at most octupoles. As explained at the end of Section 2, the B = 7 Skyrme soliton 
is believed to have octupoles as leading multipoles, but there is no numerical evidence for 
leading multipoles of higher order. Unfortunately, it seems very difficult to rule out this 
possibility in general. 

Even if one could prove (or circumvent) the assumption concerning multipoles, the exis- 
tence of attractive forces between Skyrme solitons is not sufficient to establish the inequal- 
ity C PL .19Q for infima. Physically it seems reasonable that the existence of attractive forces 
should imply the existence of minima in every sector. However, we have not been able to 
develop this observation into a mathematical proof. Further thoughts and speculations in 



this direction can be found in |T7 
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A Spherical Harmonics 

In this appendix we derive the relation between the standard spherical harmonics and the 
following functions on IR 3 — {0} used in the multipole expansion in Section 4.1 

9"dt"i}) if n>0 

Here N >0 and —N < n < N and d = \{pi — id 2 ). These functions are harmonic in their 
domain: 

AF Nn = 0. (A.2) 
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They are also homogeneous of degree — (N + 1) so that they can be written as 

1 



Nn 



r N+l 



®Nn(0,<f), (A.3) 



where (9, <p) are the usual spherical coordinates on the two-sphere centred at the origin. 
Since the Laplace operator takes the following form in spherical coordinates 

A 1 9 2 1 A ,. A . 

A = -—r + -A UJ , (A.4) 

where A u is the Laplace operator on the 2-sphere of unit radius, it follows from (|A.2|) that 

A^ Nn = -N(N+l)$ Nn . (A.5) 

Define the generator of rotations about the 3-axis 

Js = -i£ (A.6) 

and express it in terms of complex coordinates z = X\ + ixi and complex derivatives in the 
X\X<i plane: 

J 3 = zd- zd. (A.7) 

The operator d acts as a raising operator and the operator d acts as a lowering operator 
for J 3 : 

[J 3 , 8} =5 and [J 3 , 8} = -d. (A.8) 

Thus if <f) n is a function on IR 3 which is an eigenfunction of J3 with eigenvalue n then 
d<f) n is an eigenfunction with eigenvalue n + 1 provided it is not zero. Similarly d<p n is an 
eigenfunction of J3 with eigenvalue n — 1 provided it is not zero. 

It follows from (|A.l ) that 



d n F N _ nfi if n>0 
<9- n Fjv +rii0 if n < 0. 



(A.9) 



Noting that, by rotational symmetry about the 3-axis, 

J 3 $ N0 = Q for all N, (A. 10) 

we conclude that 

J 3 F Nn = nF Nn . (A. 11) 

It follows that 

® Nn = r N+1 F Nn (A.12) 
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also satisfies 

J 3 <5> Nn = n<5> Nn . (A. 13) 

Thus, to sum up, the $A? n are functions on S 2 , which are eigenfunctions of both the 
Laplace operator ( |A.5| ) and the operator J3 with eigenvalues respectively —N(N + 1) and 
n. It follows from standard harmonic analysis on S 2 that they must be proportional to the 
spherical harmonics Y^n- 

In the case n = we can determine the proportionality constant by evaluating both $jvo 
and Yn on the positive 3-axis. With the usual normalisation |jTB[ we find 



r2N + T(-i) N 

Y »° = v J n\ (A - 14) 

The relation between Y^ n and $ Nn for n 7^ is harder to compute. Let us assume initially 
that n > 0. Starting with the standard expression for the associated Legendre function in 
terms of Legendre polynomials 

P N (cosO) (A.15) 

a cos a J 

and the expression of the spherical harmonic in terms of associated Legendre function (see 
e.g. [IIH] p. 99) we have the relation 



YM = (-iryj^j fess) < A - 16 > 

Then using ( |A.14j ) and the definition of $tvo we deduce 



y Nn(e , v) = r^^tl^J^,^ ( I ) , (A , 7 ) 



4tt iV! V (JV + n)! y rf V r 
where 

De= l ' ^ +d3 (A.18) 

r a cos u pop 

and 2 = pe %ip with p = r sin # as in the main text of the paper. Then we use the commutation 
relation 

pop 

to move Dq past 83 in (|A.17 ). Noting that Dqt = we find 



w ^ = ^y^±i „(-'>:: >« ri , . (A ,o) 



4vr v /(iV-n)!(A^ + n)! VP#P 
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Now exploit that on any function / which only depends on p and x 3 



~M = -df(p,x 3 ) (A.21) 
p op z 



to conclude 



J 

Thus we finally arrive at promised relationship between Y Nn and &Nn, valid for n > 0: 



2N+1 -1 ) N + n 2 n IK x 

y^ n = J— ^ ^ ^=$jVn- A.23 

V 4tt v /(iV-n)!(iV + n)! 1 ' 

To deduce the corresponding result for n < we note that $at„ = $jv(-n) and Y/vn = 
(-l) n FiV(-r t ) for n < 0. Thus for n < 0: 



v Nn = \ t y l ) z A _24 

1 4tt v /(iV-n)!(iV + n)! V ; 
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